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We discuss the influence of atomic thermal motion on the efficiency of multimode quantum mem¬ 
ory in two configurations: over the free expand of atoms cooled beforehand in a magneto-optical 
trap, and over complete mixing of atoms in a closed cell at room temperature. We consider the 
high-speed quantum memory, and assume that writing and retrieval are short enough, and the dis¬ 
placements of atoms during these stages are negligibly small. At the same time we take in account 
thermal motion during the storage time, which, as well known, must be much longer than durations 
of all the other memory processes for successful application of memory cell in communication and 
computation. We will analyze this influence in terms of eigenmodes of the full memory cycle and 
show that distortion of the eigenmodes, caused by thermal motion, leads to the efficiency reduction. 
We will demonstrate, that in the multimode memory this interconnection has complicated character. 

PACS numbers: 42.50.Gy, 42.50.Ct, 32.80.Qk, 03.67.-a 


I. INTRODUCTION 

Over the last decade various protocols of quantum memory, based on the interaction between signal and driving 
light pulses with an ensemble of immobilized atoms were proposed [IHi. Certainly, the approximation of motionless 
atoms is natural, when one implements a quantum memory protocol on the impurities in crystals [^. There atomic 
motion is restricted by the nodes of a crystal lattice, and any spatial fluctuations are negligibly small. Besides, this 
approximation provides us with reliable results, if a full memory cycle from the beginning of a writing stage and 
until the end of a signal retrieval (including the storage time) is short enough, and the root-mean-square velocity is 
relatively small Q- 

The approximation of motionless atoms plays a significant role in analysis of multimode memory process because it 
allows us to follow the conversion of a time dependance of signal into a spatial distribution of a collective spin at the 
writing stage. Then during the storage time this spatial distribution keeps unchangeable, and after that, is converted 
back into a time profile of signal at retrieval. 

However, since the purpose of quantum memory is long-term storage of information in the generated spatial coher¬ 
ence mode, then to adequate assessment of the potential of such protocols we have to introduce thermal motion of 
atoms in a model and estimate the influence of ’’blur” on a spatial distribution of collective coherence. 

The experiments, related with light slowing, storage and manipulation on the cells with warm atomic vapors IMl, 
are quite attractive for their comparatively simple treatment. Notably, it is much easier to create an atomic ensemble 
with big number of atoms, when there is no need in its deep cooling. Moreover, a concentration of atoms in the 
ensemble can be well controlled by adjusting its temperature [T5 |. An important characteristic of the memory is its 
scalability [a [ 11 . In this respect, the ensembles of room temperature atoms are more promising than cold ones, 
which require extended cooling apparatus. 

Thermal motion of atoms is reflected in two phenomena. First of all, it is the Doppler light shift. As a result, the 
medium has a Voigt absorption profile (the convolution of Gaussian and Lorentzian profiles), what can dramatically 
change the character of the light-matter interaction due to its inhomogeneous nature. However, because of collective 
properties of the atomic ensemble, the Lorentzian profile would be determined not by the spontaneous decay rate of a 
single atom 7 , but by the decay rate of the ensemble as a whole dq, where d is the optical depth. It is the width of the 
profile compared with the Doppler width. That is why in the memory process the influence of thermal motion on the 
absorption profile would be significantly suppressed. Furthermore, the broadening of the two-photon spin transition 
can be eliminated by working in collinear geometry (e. g., co-propagating signal and driving fields) because here the 
transition is caused by two consequent processes - the absorption with following the re-emission, which have equal 
but opposite frequency shifts. 

The second phenomena induced by thermal motion, is a dependence of time of the collective spin. It occurs, because 
replacement of atoms leads to reshaping of the spin distribution, which can significantly vary the memory process 
on each stage ~ writing, retrieval, and especially storage. In particular, thermal motion would change the spatial 
coherence modes at storage, that would be reflected on an optimization mechanism of the memory efficiency | 16 l - l 20 l| . 

In this article we will study the influence of atomic motion on the efficiency of multimode quantum memory. We 
will assume that writing and retrieval are short enough, and therefore we can neglect any displacements of atoms 


2 


on these stages. At the same time we will take in account thermal motion during the storage time, which, as well 
known, must be much longer than durations of all the other memory processes for successful application of memory 
cell in communication and computation. We will analyze this influence in terms of eigenmodes of the full memory 
cycle and show that distortion of eigenmodes, caused by thermal motion, leads to the efficiency reduction. We will 
demonstrate, that in the multimode memory this interconnection has complicated character, and prove it with the 
numerical calculation. 

We will discuss two configurations of memory cells with atomic motion, related with the different experimental 
Mproaches. The first one assumes, that a signal is mapped on cold atoms, prepared in a magneto-optical trap 
. At the beginning of writing the magneto-optical trap is switched off, and the atomic cloud begins to 
expand freely. The mean temperature of the ensemble is much higher than the degeneracy temperature, so that the 
atomic motion satisfy Maxwell-Boltzmann statistics [1^. The interaction between atoms as well as their collisions 
with the walls of the cell can be neglected in the assumption that the mean free path of each atom in the cloud is much 
higher than its average displacement during the storage time Tg. Such situation occurs, for example, in experiments 
[2(1 l27j |. where authors explore the quantum memory protocol on cesium vapors with the concentration about 10® 
particles in mm^ and the average temperature about 100 . 

In the second configuration we consider a room temperature atomic vapor inside the cell extended in the longitudinal 
direction (along the signal and driving pulse propagation direction) and narrow in the transverse direction, so that 
transverse degrees of freedom are absent. We assume that during the storage time all atoms in the cell are mixed 
completely, so that ’’which atom” information completely erases and the spatial spin distribution formed on the 
writing stage becomes uniform. Such configuration of cells (with spin preserving coating deposited on the walls) 
was proposed in [l^ . where authors investigate experimentally the Raman-type quantum memory based on room 
temperature Cs-atoms. Note, the authors of cited article placed the cells inside the cavity that eliminated the spatial 
aspect there. 

The article is organized as follows. In Section [IT] we discuss the physical model of high-speed quantum memory 
and derive the main equations and their solutions for writing, storage and readout stages. Then, in Section ITIt l 
the eigenfunctions of full memory cycle are analyzed as well as the corresponding eigenvalues. Here we define the 
response functions, which describe a spatial distribution of coherence formed at the end of the write process, when 
one of the eigenfunctions is incident on the cell input. Section HIlC is devoted to the study of distortion of response 
functions during the storage due to the thermal motion of atoms. Following response functions allows us to estimate 
numerically the mobility of atoms calculating the overlap integrals of excited spin modes at the beginning and at 
the end of storage. In Section Hill we consider another configuration of the experiment associated with the mixing of 
atoms in the cell at room temperature. In Section EYl we optimize the full cycle of the memory taking into account 
the thermal motion of atoms at the storage stage. 


II. QUANTUM MEMORY WITH SLOW EXPANSION OF ATOMS 

A. Model 

Let us begin our research with the quantum memory protocol, based on the resonant interaction of a homogeneous 
ensemble of three-level atoms in the A-configuration with signal Eg and driving Ed light pulses. Atoms are situated 
inside a plane layer with length L, which is orthogonal to z-axis. We consider signal and driving light fields as plane 
waves propagating parallel to z-axis. We assume that the driving pulse is a strong classical field, and the signal pulse 
is a weak quantum field. 

Fig[T] shows stages of the full memory cycle: writing at 0 < t < T^, , storage at t < {T^, -\- Tg) and retrieval at 
(Tu, +Tg) <t < (Ttju -1- Tg -|- Tr), where T^,, Tg are durations of the each stage, correspondingly. 

We choose the durations of the signal and driving pulses much shorter than the time of spontaneous relaxation of 
the excited state. In particular, such assumption allows us to take into account the spontaneous decay only during 
the storage. This choice of interaction times corresponds to so called high-speed memory regime [^. 1^. 

Before the memory process, the ensemble of atoms is prepared in the ground state |1). During the writing stage 
in two-photon interaction the weak signal field transfer a part of atoms from the state |1) to the excited state |3), 
and at the same time the strong driving field, acting simultaneously on the supplementary transition, transfer these 
atoms from the state |3) to the ground state |2). Thus, the coherence between levels |1) and |2), which carries all 
quantum-statistical properties of the signal light pulse, is built. 

The ideal storage implies, that the coherence between levels |I) and |2) remains unchanged. However, we will 
consider two factors of its distortion: thermal motion of atoms, and spontaneous decay of the residual population 
from the excited level |3) to the ground level |1). 
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Figure 1: Sketch of the full cycle of high-speed quantum memory: a) writing, b) storage with thermal motion, c) retrieval, d) 
atomic energy levels with signal Es and driving Ed light fields. 


At the retrieval stage atoms from the level | 2 ) under the action of the driving pulse return to the level |1) through 
the excited level |3). As a result the emission of an output signal occurs, and this pulse completely reproduce the 
quantum state of the initial signal or, in the imperfect case, carry some of its properties. 

We suppose the interaction processes between the atomic ensemble and the light pulses are short enough, so 
displacement of atoms during these times are negligible and we shall consider an atomic thermal motion only during 
the storage. 

In the dipole approximation the light-matter interaction on writing and readout stages is determined by the following 
Hamiltonian: 

V = -'^dj{t)E{zj,t), (1) 
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E{zj,t) = Esizj,t) + Ed{zj,t), 

where dj is the electric dipole operator of the jth atom, located in Zj at t = 0 . 

The signal and driving fields are given as follows (in the plane wave approximation); 


Es{z,t) = + +Ac. , 

( 2 ) 

Ed{z,t)= + + 

(3) 


where ks and kd are wavenumbers of the signal and driving fields. The amplitudes d{z,t) of the signal field and a of 
the driving field are normalized so that {d^{z,t)d{z,t)) and |ap give the mean photon fluxes in photons per second 
for the light beam of area S. The operators d{z,t) and d^{z,t) obey the commutation relations; 

[d{z,t),d\z,t')] =S{t-t'), (4) 

[d{z,t),d\z',t)] = c(^l-Hz-z'). (5) 

It is natural to describe the light-matter interaction in terms of collective operators, which can be defined as a 
superposition of the microscopic variables: 

dmn{z,t) ='^\m)(n\iS {z - Zi), (6) 

i 

Nmiz,t) ='^\m){m\iS {z - Zi), (7) 

i 

= (^Nm{z, t) - Nn{z, t)^ d(z- z') . 




( 8 ) 
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Such description allows us to use a continuous variable associated with the observation point instead of discrete 
coordinates of atom positions. Although the collective operators turns out to be sufficiently singular, it was shown in 
[ 2 ^ how to perform its smoothing. 

Let us rewrite the Hamiltonian in the collective variables: 


V= dz 


ih[gd{z,t)a- 3 i{z,t) — h.c.) + ih{p,ds 2 {z,t) — h.c.) 


ikdZ 


(9) 


Coupling between the light pulses and the atomic ensemble is defined by the coupling constant g and the Rabi 
frequency which we choose real for simplicity: 

1/2 / ^ 1/2 


g = 


2eoficS 


^31) 


n = 




2€oficS 


l32- 


( 10 ) 


Here d^i and d ^2 are the respective matrix elements of the electric dipole operator. 

It is important to notice, that the probability of absorption of signal photons on the transition |1) — |3) is determined 
not by its spectral width 7 , but by the product of 7 and the optical depth d. That is why, in spite of the small 
probability of absorption by a single atom, the probability of absorption by the whole ensemble is high. 

Since the mean number of photons in the signal pulse is small, we consider the population of the level |1) constant 
during the full memory cycle and treat it as a real number N [281 l29j|. Taking this into account, it is convenient to 
re-normalize the coherences and introduce new atomic operators b{z,t) and c{z,t)\ 


b{z,t) = ai 2 {z,t)/VN, 
c{z,t) = ai 3 {z,t)/y/N, 


( 11 ) 

( 12 ) 


which satisfy the bosonic commutation relations. For further calculations we shall derive the simultaneous commuta¬ 
tors: 


b{z,t),b'^iz',t) =S{z-z'), 
[c{z,t),c'^{z',t)] = 5{z - z') 


(13) 

(14) 


Rewriting the Hamiltonian with the new operators, we define a new coupling constant g^, which is the product of g 
and '/N: g^ = '/Ng- This form of the coupling constant emphasize, that the light-matter interaction is determined 
by the number of atoms in the ensemble, and collective effects become apparent for large N resulting in strong 
interaction. 

Let us derive a closed set of Heisenberg equations for operators d[z, t), b{z, t), and c(z, t) that describes the evolution 
on the writing and the retrieval stages, i.e. at 0 < t and (T^, -|- Tg) <t< {T^ -I- Ts -|- Tr): 

d 

-d{z,t) =-gNc{z,t), (15) 


dz 

m 


c{z,t) = gN d{z,t) + nb{z,t), 
b{z,t) =—ftc{z,t). 


(16) 

(17) 


The derivation of the equations and the way of their solution can be find in [28|, l29| , where the protocol of high-speed 
quantum memory on cold atoms were proposed. Solutions (in dimensionless variables) connect the amplitude of the 
signal field dmit) on the input of the cell with the coherence b(z,T^) at the end of writing. 


bi^z^T^') — / dt Gab(^z ^ unc, 


(18) 


Gab{zj) = -^J^ di'gab{z,i')gl^,[z,i-i'), gab{z,i) = e **Jo 0u,(t), 

as well as the coherence b{z,Tw + Tg) at the end of the storage with the output field dout(t) obtained in backward 
retrieval 


doutik) — / dz Gbai^z^ t')b{^z^ Tyj Tg) -t- vclc^ 

Jo 

1 /■* _ . . . 
Gba{z,t) = dt'gba{z,t')gl^{z,t-t'), 


( 19 ) 


gba{z,i) = e 0r(t), 
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Here Jq is the zero order Bessel function of the first kind, Qw{i), Qr{i) are the window-functions (they are equal 
to one during the writing and the retrieval stages, correspondingly, and zero-valued for all the other time intervals). 
Dimensionless coordinate z and time t defined as follows: 

z = t = tn, (20) 

Henceforth we will omit ’’tilde” over variables, regarding them dimensionless, if the otherwise not stated. 

Note, that in Eqs. (Ulii-dlll) there are terms, which marked as vac and correspond to contributions of the subsystems 
in the vacuum state. There is no need to specify them, since their contributions vanish when normally ordered 
expectation values are taken. 

To consider atomic motion during the storage and estimate its influence on the stored coherence b(z,T^), let 
us introduce a subensemble of atoms moving co-directionally with a velocity Vz- Then evolution of the coherence 
b{z,Tuj',Vz) of this subensemble during the storage is determined by the following equation 

The transition from the ensemble to the subensemble could be formally done as far back as in the Hamiltonian o, 
if one redefine properly the collective variables b{z,t) and c{z,t) for the subensemble. However, we prefer to perform 
it now to underline, that we consider atomic motion only during the storage. Note, the equations (nsi-dni) and their 
solutions dll-dill) are valid not only for the ensemble of motionless atoms, but also for any its subensemble, which 
moves as a whole with some certain velocity. 

The solution of Eq. m at the end of the storage -|- Eg can be got in a simple form: 

b{z, Tzu + Ts; Vz) = b{z - VzT^, T„; u^), (22) 

where Vz is a dimensionless velocity defined according (1^ . The equation (1^ shows that at the end of the storage 
stage the domain of function b{z, + Ts]Vz) is shifted on VzTg, from z € [0, L] at the beginning of the memory cycle 
to z € [t'zEg, L + VzTs] at the end of the storage. At the same time, the profile of the function remains unchangeable. 
Thereby, to describe the motion of all atoms, we shall take into account all contributions b(z,Ty, + Ts]Vz) from each 
subensemble to the coherence b{z,Tyj + Tg) of the whole ensemble, and this coherence is the initial condition for the 
retrieval stage. 

Further analysis of the solutions derived in this Section we will carry out in terms of eigenfunctions and eigenvalues 
of the memory cycle. 


B. Eigenfunctions of the memory cycle and response functions 

Let us start with the consideration of cold atoms, which stay immobilized during the full memory cycle. The 
connection between the input and the output (retrieved) fields is described by the integral transformation: 

Tu, 

aout{t) = J dt' ain{T^u-t')G{t,t')-\-vac. (23) 

0 

Here G{t, t') is the real kernel of the integral operator. It can be expressed by the kernel of the writing stage G'at,(z, t), 
converted the input field ain(t) to the spin coherence b{z, Ty,), and the kernel of the retrieval stage Gba(.z, t), converted 
the spin coherence b{z,Ty,) to the output field aout(t): 

L 

Git, t') = J dz Gabiz, t)Gba{z, t'). (24) 

0 

Since in the high-speed memory regime Gabiz,t) = Gbaiz,t), than the kernel G(t,t') is symmetric with respect to 
permutation of the arguments t t'. This means we have a right to derive the equation for its eigenfunctions 

and eigenvalues in the form: 

\/A'<(>j(t) = j dt' G{t,t')(j)i{t'). 

0 


(25) 
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Figure 2: The first five eigenvalues as a function of the duration T of the writing and the retrieval processes, when Tr = = T 

and L = 10. The red columns correspond to the eigenfunctions with efficiency r]i > 50%, the blue columns - with efficiency 
Vi < 50%. 


The functions form the complete orthonormal set: 




dt ; 

0 

(26) 

OO 

2 = 1 

(27) 


One should note, that such definition of the eigenfunctions assumes the equality of the writing and the retrieval 
durations. In more general case, the duration of the retrieval can exceed the writing time, so the arguments of the 
kernel will be defined on the different domains. However, even in this case it is possible to symmetrize the 

kernel and find its eigenfunctions and eigenvalues (see Appendix A). 

One can derive the equation equivalent to the formula (I25L representing the kernel G(t,t') as a bilinear quadratic 
form of its eigenfunctions, where the corresponding eigenvalues are decomposition coefficients: 

OO 

= (28) 
i=l 

Such representation commonly called the Schmidt decomposition, and 4>i are referred to as the Schmidt modes. 

The diagram in Fig. [5] depicts the first five eigenvalues as a function of the duration of the writing and the retrieval 
light pulses, when Tr = = T and the length of the atomic layer is L = 10. It shows the eigenfunctions with 

the efficiency higher then 50% (red columns) that means quantum memory regime, and below it (blue columns). In 
other words, if the input signal profile coincides with the zth eigenfunction, than the quantum memory performs like 
a beamsplitter with the transmission coefficient rji, which equals to the squared ith eigenvalue and determines the 
efficiency of the memory: 


Jo dt{a\^^{t)aout{t)) , , . 

m = -= Xi. (29) 

fo dt(aj^(t)ain(t)} 

Let us make a remark, that if a time profile of the input pulse is not identical to any of eigenfunction, then the 
description of the protocol as a beamsplitter with a certain coefficient is not valid longer. The properties of such 
protocol will be determined by the set of numbers associated with projections of the input temporary profile on the 
eigenfunctions of the memory [s^. 

As one can see from the picture, for any T eigenvalues decrease rapidly and, in fact, only the first two are noticeably 
different from zero. That is why we will consider further only the first two eigenfunctions. 

Note, this memory model is valid only when T < L [^ . 
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Figure 3: The eigenfunctions (jiiit) and < 7 ^ 2 (t) of the memory cycle for motionless atoms (left column, frames a and b), and their 
squares (right column, frames c and d). Dimensionless time t is given in units of 


Further we will consider only the case of backward retrieval at L = 10 and T = 5.5. We exploited such values of 
parameters in [30l | to find eigenvalues and eigenfunctions in the case of motionless atoms, and now we want to compare 
them with the case of thermal motion. 

Fig. [ 3 ] shows the first two eigenfunctions of the full memory cycle for motionless atoms (left column) and their 
squares (right column). 

Let us now consider a situation, when the signal with temporary profile of one of the memory eigenfunctions 
incidents on the input face of the atomic cell, and find out, what kind of ’’response” of the medium would caused by 
such a field. We shall call this transformation the ’’half-cycle” as the opposite to the full memory cycle of writing and 
retrieval of the signal: 


y/JTiiJiiz) = J dtGab{t,z)(l)i{t). (30) 

0 

As will be proved below, the set is orthonormal and the normalization factor ^/Jli is related with the eigenvalues 
Ai of the full cycle. We shall call products the response functions of the medium. 

In the contrast to the kernel of the full memory cycle, the kernel Gab{t,z) is not symmetric in respect to 

the permutation of its temporary argument t and spatial argument z, and so the Schmidt decomposition is ineligible 
here. However, due to the completeness of the eigenfunctions {(/>i(t)}“i, we can expand Gab{t,z) into the series: 

00 

Gabit, z) = '^ gi{z)(l),{t), (31) 

i=l 

here gi(z) are the expansion coefficients. Now we can find the relation between gi{z) and y/JIlipi{z). For this we 
multiply the right and the left parts of Eq. (1311) on (j)j(t) and then integrate this expression by time from 0 to T^,. 
Taking into account the orthonormality of {4'i{t)}^i, we get 

= 9j{z)- (32) 

This relation means, that the ith response function of the medium yfjlfifii^z) is the ith expansion coefficient of the 
Gab{t, z) by the eigenfunctions , so the equation (I5T|) can be rewritten in form 

00 

Gab{t, - 2 ) = X! (33) 

2 = 1 


Now let us prove, that functions 'il^iiz) compose a complete orthonormal set, and find the normalization factors 
































































Figure 4: Squares of the response functions. Top row - | AL)p, bottom row - | AZ/)p. In frames a) and 

b) AL = 0 (motionless atoms), in frames c) and d) - AL = 2, in frames e) and f) - AL = 10. Dashed curves are given for 
ease of comparison and correspond to AL = 0. Dimensionless coordinate 2 : is given in units of the effective optical depth (see 
subsection III C 2ll . 


Express the scalar product of the zth and the jth response functions, using Eq. dsni): 

pL pL pT-u, p'T'^ 

dzil}i{z)il}j{z) = / / / dzdtdt'Gab{t',z)Gab{t,z)(l)i{t')(j)j{t) (34) 

Jo Jo Jo Jo 

Taking into account the Eqs. we obtain 

/■L 

^/JTiy/JTj / (35) 

Jo 

This means that functions orthonormal, with the normalization factors = jAi. As a result the following 

expansion for Gab(t, z) can be written: 


Gab{t, -z) = E (36) 

i=l 

The completeness of the set {ipi{z)}^i is followed from the linearity of the integral transform (I30|) and the completeness 
of the set {^!'i(t)}“i- 

The functions ’4>i{z) can be conditionally called the eigenmodes of the spin system. Further, we will be mostly 
interested not in these eigenmodes, but in the corresponding response functions •^^'ipi{z), squares of which reveal 
the distribution of the excitations within the medium. 

Figure |4^ shows the spatial dependence of the first and second response functions for the motionless atoms. The 
area under the curve for the first response function is ^/Xl = 1.0, for the second one - = 0.8 that is well agreed 

with the corresponding storage efficiencies. Besides, one can take notice, that the first curve has only one peak at 
2 ; = 0, at the same time the second curve has two peaks at 2 ; = 0 and at z ~ 3.9. Thus, in the first case the excitations 
are localized mainly near the input face of the cell, whereas in the second case they partially ’’pushed” in the middle. 
In the next Section we will follow how the obtained curves changes due to the thermal motion at the storage stage. 


C. Thermal motion at the storage stage for slow atoms 

1. The limits of applicability 

Now let us consider the atomic motion in the longitudinal direction at the storage stage, when < t < + Tg. 

Well known, that an ideal gas can be described in the framework of classical statistics, when its temperature exceeds 
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the degeneracy temperature [25j |: 


r» 


Smk ’ 


(37) 


here T is a temperature, n is a volume concentration of the atoms, to is a mass of a single particle, k - Boltzmann 
constant, h - Plank constant. The macroscopic behavior of nondegenerate gas obeys Maxwell-Boltzmann statistics. 

We consider the case, when an atomic ensemble has Maxwell speed distribution, because it covers wide class of the 
quantum memory systems. For example, it takes place in the works [^. where authors experimentally consider 
quantum memory on the caesium vapors with the average concentration of atoms about 10® particles in 1 mm® at 
the temperature about 100 fiK. One can verify, that such parameters satisfy the inequality (1371) . 

In Section III Al we derived the equation (I22L which bonds a coherence of spin subensemble moving as a whole in 
the longitudinal direction with the velocity Vz at the moment t = T^j and a coherence of the same subensemble at the 
end of the storage t = + Tg. Let us rewrite this expression for the ith response function: 




Vz ,t —Tiu+Tg 


= - VzTs) 


Vz ,t—T.w 


and pass from the subensemble to the whole ensemble with the Maxwell speed distribution: 


\/%ipi{z) 




^ + oo 


dvz e - VzTg) 


Vzft — Tyj 


(38) 


(39) 


where Uz is the root-mean-square velocity of the particles in the longitudinal direction. The mean extension of the 
spin ensemble during the storage is given by AL: 


AL = TgUz, (40) 

which specify the mean temperature of the ensemble (for the given Tg). Hereinafter we will use denotation AL) 
instead of ipi{z) 


t=T„+T, 


2. Scaling of the coordinate 

One can see that the definition of dimensionless coordinate (EUl) for motionless atoms implies the atomic concen¬ 
tration N as constant, then the dimensionless coordinate coincides with an effective optical depth (which differ from 
the real optical depth in H/y times). In other words, in the case of motionless atoms the dimensionless coordinate 
is measured in units of the effective optical depth. When the dimensional coordinate varies in z S [0,L], the dimen¬ 
sionless coordinate varies in z G [0,L]. (In this subsection we again return to the old notations with tilde for the 
dimensionless variables.) 

As a result of thermal motion, the homogeneous distribution of atoms is disturbed, and the concentration of atoms 
N becomes a function of z, with a domain z G [—oo,-|-oo] (nevertheless, the most of atoms are located within 
z G [—AL,L -|- AL])). This change breaks the direct correspondence between the dimensionless coordinate z and 
optical depth. We want to redefine the dimensionless coordinate so that to restore this equivalence. We will introduce 
new dimensionless coordinate z, which will vary in the same interval z G [0,L], as the old one. By this reason we 
carry out the scaling of the coordinate. 

The main idea of the scaling is to find the relation z = /(z) between the old dimensionless coordinate z and a new 
one z such that concentration N[z) is again homogeneous. 

As far as the longitudinal motion does not change the whole number of atoms in the system, we can always write 
the following equality 


dzN{2)= / dzN{z). 


(41) 


Here the left part corresponds to the whole number of atoms before storage, and N{z) is a constant. The right 
part is the same number after the thermal motion, N can be obtained from N{z) by averaging it over the Maxwell 
distribution, and we demand N{z) is also constant. Let us replace in the right part z on z: z = /(z), dz = f'{z)dz, 
and L = L, we get 


‘ dzN{i)= r dz f{~z)N{f{i)). 
Jo 


(42) 
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From here one can obtain 


f'Cz) 


Nji) 


(43) 


As far as N{z) and A^(/(5)) are known, we can solve this equation numerically, and reconstruct the scaling function 

m- 


3. Response functions for mobile atoms 

Thermal motion of atoms leads to the spatial redistribution of the stored spin coherence. Let us now apply the 
scaling procedure described in the previous subsection, and follow the influence of such redistribution on the response 
functions AL). Hereinafter we again omit ’’bar” over z, regarding it as the dimensionless variable expressed 

in optical depth units (with recipe given above). 

Figs. IUd and 131: show the spatial dependence of squares of the first and the second response functions for mobile 
atoms, when the mean shift of atoms equals AL = 2 (i.e. one fifth of the cell length) and AL = 10 (i.e. the full cell 
length), correspondingly. As before, the top row (blue curves) corresponds to the first response function, the bottom 
row (purple curves) - to the second. For simplicity of comparison, the dashed curves mark the case of motionless 
atoms. First of all we see, that for small displacements {AL = 2) the shapes of the curves and the areas under the 
curves have changed slightly comparing with the case of motionless atoms. So the spatial distribution of excitation in 
the medium changes insignihcantly, and we can expect that after retrieval we receive an output held with a time prohle 
similar to the corresponding input eigenmode with high efficiency. Note, the both curves have the same number of 
peaks, as they had in the case of motionless atoms. However, due to thermal motion these peaks noticeably subsided, 
and the peak at z ~ 3.9 for the second curve shifted to the right in z r: 4.6. From the hgure one can conclude that the 
stored excitation in the medium was ’’blurred” by the thermal motion over the whole spin ensemble. This is clearly 
seen in the case of the mean displacement equals AL = 10. Besides, it shows, that the higher number of zeros of the 
mode 'i/’i(^) and the lower ’’contrast” of its peaks, the faster ’’blur” will destroy this mode. In particular, we see, that 
at AL = 10 the peak at the beginning of the first response function, in spite of its subsidence, is preserved, whereas 
the both peaks of the second response function have vanished. 


4- Overlap integrals 

We have been convinced, that thermal motion redistribute the stored spin excitations from the initially excited mode 
over the others spin modes. At the same time, the main role in the recovery of the signal after such redistribution will 
play the modes with the highest eigenvalues y/Xi. Therefore one need to retrieve the signal not only from the initially 
excited mode but also from all the others with high y/Xi for the effective retrieval. To confirm such interpretation 
quantitatively, we introduce the overlap integrals Qij{AL) of the ith ’’blurred” mode of the thermal atomic ensemble 
with the jth eigenmode of the immobilized ensemble: 


Qij{AL) = [ dz^i{z-AL)ij;j{z;0), (44) 

Jo 


where ipj{z] 0) = ipj{z)- It is seen from the definition that overlap integrals can be rewritten in the matrix form, which 
has Qij as the element of the ith row and the jth column. 

Then, for AL = 2 the overlap integrals are given by 


Q 


fO.92 0.1l\ 
1^0.11 0.74:) ’ 


and for AL = 10 we got 


/0.68 0.32\ 
1^0.32 0.39j ■ 


These numerical calculations confirm our previous conclusions. The square of the first response function has only 
one peak, whereas the square of the second one - two small peaks. As a result, the second response function will be 
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changed by thermal motion more significantly than the first one, what can be seen from comparison of Qn and (522- 
In other words, the first function can be stored better, than the second. The equal values of Q 12 and Q 21 reflects the 
number of excitations, which flowed from the first mode to the second one and vice versa. One can notice, that in the 
case of AT = 0 (motionless atoms) Q 12 and Q 21 are equal to zero. Thus, in any cases the matrix Q is symmetrical. 

One can evaluate the efficiency of the full memory cycle over the overlap integrals. For this purpose let us expand 
the ’’blurred” response function AL) by the full orthonormal set Then, taking into account 

Eq. dSl) we get 

(30 

y/xltljiiz-, AL) = (45) 

1=1 

Using this formula, Eq. m, and Eq. dMl), we can derive the time profile (t) of the output signal at the readout, 
when the incident light pulse had the profile of the *th eigenfunction: 

00 

= (46) 

1=1 

Substituting Eq. (l46ll into the definition of the memory efficiency (l29ll . we got 

(XD 

(47) 

1=1 

Using this formula, we obtain the following values of efficiency. In the case of the mean displacement AL = 2 for the 
input pulse shaped like the efficiency of the full memory cycle is r]i = 87%, like (j>2{t) - the efficiency 772 = 46%. 

In the case of AL = 10 we obtained ryi = 55% and 772 = 22%. All these values completely coincide with the results of 
direct numerical calculations of 77^, based on the initial integral transformations. As one can see, even for relatively 
large temperatures of atoms and without any additional optimizations, one can choose the time profiles, which provide 
the quantum level of the efficiency (higher than 50% ). 

Note, the diagonal elements of Q correspond to single-mode efficiency (so called beamsplitter-efficiency) of quantum 
memory, when one readout the signal from the same eigenmode, in which it was written. With respect to this mode, 
quantum memory works as a beamsplitter with transmission coefficient equals to storage efficiency. Single-mode 
efficiency is remarkable because it describes not only the quality of storage of the photon numbers, but also all the 
other moments of quantum state distribution of the field. It is important to remember that this is the only case where 
the efficiency is a universal characteristic of the system. If the recorded signal is not an eigenfunction of memory, or 
reading mode does not match the writing mode, the relationship between efhciency and preservation of other quantum 
properties of the field becomes more complicated [3l| . 

In Fig. Othe squares of are plotted. As before, the top row (blue curves) corresponds to the first eigenfunc¬ 

tion, the bottom row (purple curves) - to the second. Dashed curves correspond to AL = 0 and given for comparison. 
It is interesting to notice that the redistribution of excitation over the full memory cycle eigenfunctions looks more 
visual, comparing with the same process for the response functions (see Fig. |3|). In particular, one can see the second 
peak, specific for which grows on |(('i“*(t)P at AL = 10. 


III. STORAGE IN THE ATOMIC CELL AT ROOM TEMPERATURE 

Let us now consider the atomic ensemble at room temperature. In this case we can not longer discuss mean 
displacements of atoms of the order of cell length, as it was in the previous Section. We will suppose that at room 
temperature the full mixing of the atoms in the cell during the storage time occurs, and at the end of the storage the 
distribution of spin coherence becomes homogeneous. 

As opposite to the model discussed before, we will suppose now that the atomic ensemble is situated not in a free 
space, but in the closed prolate cell with length L (along the signal and the driving pulses propagation), and with the 
narrow cross-section (with no transverse degrees of freedom). As before, we will study the regime of the high-speed 
quantum memory and treat atoms as immobilized at the writing and readout stages, due to the shortness of the 
interaction times. 

We imply that cell coated with the spin protecting alkene coating and any atomic relaxations caused by 

collisions with walls can be eliminated. 

Because now atoms do not leave the cell, and their concentration stays uniform, the definition of the dimensionless 
coordinate remains similar to the case of motionless atoms, and we do not need to apply any scaling. 
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Figure 5: The time profiles of the intensity of the output fields. Top row - |(/)i“*(t)|^; bottom row - |(/) 2 “*(t)|^. In frames a) 
and b) AL = 2, in frames c) and d) AL = 10. Dashed curves are given for ease of comparison and correspond to AL = 0. 
Dimensionless time t is given in units of 



Figure 6: a) The square of the first response function right after the writing (dashed curve) and after the storage (solid curve), 
b) The time profile of the output field intensity stored on room temperature atoms (solid curve) and motionless atoms (dashed 
curve). Dimensionless time t is given in units of 11“^. Dimensionless coordinate 2 is given in units of effective optical depth 
(see subsection III C 21) . 


We shall emphasize, that in contrast to immobilized atoms, any spin distribution, formed in the writing process, will 
be blurred by the thermal motion and converted into uniform, independently of the shape of input pulse. However, 
the input pulse shape will determine the total number of spin excitations. As a result, the output signal shape will be 
always the same, and only difference will be revealed in the efficiency of the full memory cycle. The latter, as before, 
can be calculated with the help of the overlap integrals dUl). 

Thus, in order to analyze this problem numerically, we should examine only the retrieval, assuming that the spatial 
distribution of the spin coherence at the end of storage stage is known (this is the constant function, which level 
corresponds to the total number of spin excitations). We have confined ourselves to the case of the input held prohle 
matched with ^i(t). 

The dashed curve in the Fig. [5^ shows the square of the response function right after the writing stage. 

The solid curve on this plot corresponds to the same function, but after the storage time. In the Fig. [BJo the solid 
curve depicts the intensity prohle of the output held, and the dashed line is given for comparison with the case of 
motionless atoms. One can see, that complete mixing of atoms inside the cell leads to redistribution of the retrieved 
intensity and to formation of two peaks at t ~ 0.6 and t ~ 3.4, which are characteristic for the hrst and the second 
eigenfunctions. The hrst peak is higher than the second one, since the efficiency for the hrst mode is higher. 

The numerical calculations bring out the high values of the overlap integrals for room temperature atoms: Qn = 0.73 
and Qi 2 = 0.56 that provide the total memory efficiency ?7i=82%. Thereby, complete mixing of atoms inside the cell, 
contrary to our expectations, does not lead to the hagrant damage of the memory, comparing with the case of slow 
atoms in a free space considered before. In particular, at AL = 10 we obtained the worse results for slow atoms in a 
free space. 
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Thus, we can conclude that, for the given optical depth, the choice of a particular experimental configuration may 
lead to improvement or deterioration of the results depending on the overlap of the stored spin distribution with the 
most efficient spin eigenmodes. 


IV. OPTIMIZATION OF FULL MEMORY CYCLE 

In this article, to ensure the maximum efficiency of the memory protocol, we applied the optimization based on the 
analysis of the kernel of integral transformation and fitting the optimal signal profile. This approach was introduced 
as a tool of the quantum memory optimization in (35l |. (Note, that the other well-known optimization method based 
on driving field shaping [l^[i3|, remains outside of this discussion.) It is clear that thermal motion of atoms at the 
storage stage leads to change of the kernel of integral transformation: the eigenfunctions for the case of immobilized 
atoms will no longer be the eigenfunctions of the task with mobile atoms. Thus, it is natural to analyze the new 
integral transformation taking into account the motion of the atoms, and specify which profile of incident pulse ensures 
the maximum efficiency. 

Let us deduce a new kernel of the full memory cycle AL) instead of Eq. (IMl) . Previously, both factors 

Gab{t, z) and Gba(t, z), which describe singly the writing and readout stages, were identical; now we include into the 
first of them the atomic motion during the storage. Formally, this brings to replacement of the response functions 
by the average response function AL) into the expansion of Gab{t, z). Consequently, the new kernel for the 

writing stage Gab{t, z; AL) is given by 


OO 

Gab{t,z-,AL) = Y,^,'4,,{z-AL)4>,{t). (48) 

i^l 

The kernel of the integral transformation for the reading stage remains unchanged. As a result, the new 

kernel for the full memory cycle AL) is following: 

L 

G{t,t']AL) = j dz Gah{t,z]AL)Gha{t' (49) 
0 

Developing the kernels Gab{t, z; AL) and Gba{t',z) as series in eigenfunctions (j>i{t), one can be certain that they 
remain symmetrical with respect to permutations of the arguments t and t', and the new eigenfunction problem can 
be posted: 


_ 

^X^{AL)(j),{t-AL) = j dt' G{t,t'-'AL)(j),(t'-'AL), (50) 

0 

where AL) are eigenfunctions, and ^Jxi{AL) are eigenvalues of the memory on the thermal atomic ensemble. 
Note, these sets are different for each temperature. 

Fig. [7] shows the squares of the first (blue curve) and the second (purple curve) eigenfunctions for AL = 2 (left 
column) and AL = 10 (right column). For comparison the analogous curves for motionless atoms are plotted by 
dashed lines. It can be seen that the time profiles of intensities significantly changed. In particular, we take notice 
that profile of |<('i(t; AL)p at AL = 10 became alike the intensity profile of the output pulse for room-temperature 
atoms (see Fig. Eb). However, the main thing here is that for the new eigenfunctions the memory efficiency have 
changed considerably. Numerical calculation shows that after the optimization the new efficiencies equal rji = 94% 
and r ]2 = 41% at AL = 2, and rji = 74% and 772 = 3% at AL = 10. (Let us remind that before we got 771 = 87%, 
772 = 46% at AL = 2, and 771 = 55%, 772 = 22% at AL = 10.) Although the value 772 turns out to be lower than 
50% threshold of quantum memory, the proposed optimization procedure enabled to increase 771 significantly, even 
for large displacements AL. 


V. CONCLUSION 

We have analyzed the impact of the longitudinal thermal motion on the operation of multimode quantum memory, 
and have shown that atomic mobility yields the redistribution of spin excitations, which can be described in terms 
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Figure 7: The squares of eigenfunctions AL) and (j) 2 {t-,AL) of full memory cycle including storage stage (solid curves). 
In frames a) and b) AL = 2, in frames c) and d) AL = 10. Dashed curves are given for ease of comparison and correspond to 
AL = 0. Dimensionless time t is given in units of 


of the response functions. Redistribution of the excitations between the response functions causes the distortion of 
the retrieved field profile: it turns out to be different from the incident one, and the disparity increases with the 
temperature of the ensemble. As a result, the efficiency of the full memory cycle decreases. 

It is important to note that even with significant displacements of atoms during the storage time, the memory 
remains to be quantum. Moreover, the quantum nature of the memory is preserved even with the complete mixing of 
atoms in the cell. This once again confirms that the quantum state of light is recorded on the collective of atoms but 
not on the individual ones, and the atomic movement in this collective, although it leads to distortion of the mode 
structure of memory, but does not affect drastically the ability of the quantum storage. 

Note, in some situations the longitudinal movement can cause not a decrease of the efficiency of the quantum 
memory, but the contrary, an increase of it. For example, when one record a mode of the full memory cycle with 
a relatively small eigenvalue, thermal spin redistribution can lead to population of readout modes with the larger 
eigenvalues that should result in the increase of the efficiency. 

In this paper, we analyzed the role of the thermal motion in the model of high-speed quantum memory. The 
choice of the model is primarily caused by the ability to ignore the displacements of the atoms on the writing and 
readout stages, setting them short. However, for other models of quantum memory (adiabatic, QND, Raman), this 
statement of the problem is possible, if we assume that the atomic velocities are sufficiently small, whereas the ratio 
between the storage time and interaction time is large. Different memory models provide different mode structures. 
For example, for the adiabatic memory one can assume the minor disruption when the first mode is stored, because, 
as we have shown, the losses due to the thermal motion depend on the shape of the response function which formed 
in the medium after writing. 

We have proposed to optimize the memory taking into account the thermal motion at the storage, and fulfilled this 
optimization by finding the new field profiles, which interact optimally with a mobile atomic ensemble. 

One of the unexpected results is that the quantum memory at room temperature may be even more effective than 
with cold atoms. This is winning for the practical application of the quantum memory on the atomic ensembles, in 
particular, it is essential for its scalability. 

We can conclude that for the effective operation of quantum memory is important not only to preserve the spin 
excitations, but also to provide their optimal spatial distribution. 
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Appendix A: Eigenfunctions in the case of different durations of writing and retrieval 

According to Eq. (1^ the evolution of the signal field during its writing/readout is determined by the kernel of the 
integral equation As known, the initial set of the equations describes the system development only during 

the limited time interval, thereby we can rewrite the kernel in the following way: 

= Go(t,t')0r(i)0»(i'), (Al) 

where the factors 0r(t) and 0u)(t) are unequal to zero only during the interaction times: 

0r(t) = 1 at 0 <t <Tr and Qwit) = 1 at 0 < t < T^,. (A2) 

Here and are durations of the writing and the retrieval, respectively. The function Go{t, t') turn into the kernel 
G(t,t') at the stationary conditions, i.e. at T^,, —>■ oo. For the most models of quantum memory this function is 
symmetrical with respect to permutation of its arguments, Go(t,t') = Go{t',t). However, as one can see, for limited 
observation times the kernel G{t, t') preserve this property only at Ty, = T^. The presence of two different 0-functions 
destroys the symmetry of the kernel relative to its arguments permutation t ^ t'. Nevertheless, it is possible to 
provide the symmetry with respect to kt ^ t', where k = T^/Tr- Indeed, we can write G(t,t') = G{kt,t'), then we 
get 

G{kt, t') = Go{kt, t')Qr{kt)Q^{t') = Goikt, t')Q^,{kt)Q^{t') = G{t', kt). (A3) 

Therefore, the Schmidt decomposition for this kernel can be given as 

= G{kt,t') = Y, (A4) 

i 

where Xi and are the eigenvalues and the eigenfunctions of the kernel G, so that 

\/kXi (piikt) = f dt'(j)i{t') G{kt,t'). (A5) 

Jo 

The set of functions is determined at the domain [0,Tu,] and satisfy the conditions: 

dt(j)i{t)(j)j{t) = Sij, = 5{t - t'). (A6) 

i 

The functions = \/k 4>i(kt) are given at the domain [0,Tr], and, one can easily see, also meet the conditions of 
orthonormality and completeness: 



dtipi{t)(pj{t) = Sij, (A7) 

i 

The sets of functions {(l>i{t)}^^ and {‘Pi{t)}^i can be interpreted as eigenfunctions at the writing and the retrieval 
stages, correspondingly, and it is obvious, they can be obtained from each other by a simple scaling the time with the 
coefficient k. 
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